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Abstract 

The  intensive  use  of  the  connected  components  algorithm  in  image  analysis  and 
robot  vision  calls  for  a  very  fast  implementation  of  such  algorithm  suitable  for  real- 
time applications.  This  paper  presents  a  hardware  design  implementing  the  algo- 
rithm due  to  Schwartz,  Sharir,  and  Siegel[8].  A  prototype  board,  without  using  spe- 
cial VLSI  chips,  had  been  constructed  in  Robotics  Research  Lab  of  Courant  Insti- 
tute. It  can  compute  the  connected  components  of  a  512x512  binary  image  in  few 
video  frame  times  (about  300  ms).  A  real-time  version  (video  speed)  in  VLSI  is 
proposed. 


1.  Introduction 

The  labeling  of  connected  components  of  a  binary  image  is  a  fundamental  prob- 
lem in  image  analysis.  An  early  method  was  developed  by  Rosenfeld  and  Pfaltz  [7] 
in  1966;  it  uses  a  pair  of  arrays,  one  containing  the  current  region  label  and  the 
other  containing  its  smallest  equivalent  label.  This  algorithm  processes  the  image 
from  top  to  bottom  to  compute  label  equivalences,  storing  the  result  in  the  arrays. 
A  second  pass  reassigns  each  label  to  its  smallest  equivalent  label.  Lumia,  Shapiro, 
and  Zuniga  [3]  improved  the  previous  method  by  using  a  very  short  equivalence 
table,  which  needs  to  cover  only  one  line.  Schwartz,  Sharir,  and  Siegel  [8]  present  a 
method  which  uses  bracket  marking  to  associate  equivalent  groups.  This  method 
enables  one  to  compute  the  component  numbers  for  each  pixel  on  the  fly,  by  using 
an  relative  small  auxiliary  "bracket  table'.  More  interestingly,  this  algorithm  uses 
only  simple  form  of  pushdown-stacks,  making  high-speed  hardware  realization  pos- 
sible. In  addition  to  the  above  mentioned  sequential  algorithms,  Shiloach  and  Vish- 
kin  [9]  present  a  logarithmic-time  connected  components  algorithm  known  for  mas- 
sively parallel  computing  systems  (e.g.  1  processor  per  pixel)  connected  in  shuffle- 
exchange  or  other  similar  pattern. 

In  this  paper,  we  present  a  hardware  implementation  of  the  algorithm  due  to 
Schwartz,  Sharir,  and  Siegel[8].  In  the  following  sections,  data  structures  and  algo- 
rithm are  described  in  detail.  Then  the  hardware  design  realizing  the  algorithm  is 
presented.  A  prototype  Connected  Components  Board,  which  does  not  use  any  spe- 
cially designed  VLSI  chips,  has  been  constructed  at  Robotics  Research  Lab.  of 
Courant  Institute,  NYU.  It  can  label  the  connected  components  of  a  512x512  binary 
image  in  just  a  few  video-frame  times  (about  300  ms).  Real-time  computation  (i.e. 
video  speed)  can  be  achieved  without  significant  redesigning  by  pipelining  the  two 
major  computing  components  of  the  board.  For  a  connected  components  algorithm 
to   be  suitable   for  intensive  use  in  real-time  image  processing  and  robot  vision 


CCB  P«ge  2 

applications,  it  must  allow  for  a  very  fast  implementation.  In  the  last  section  we  dis- 
cuss this  enhancement  and  outline  a  possible  VLSI  implementation  of  the  algorithm. 

2.  The  Algorithm  and  Its  Implementation 

The  algorithm  that  we  have  implemented  was  first  presented  in  [8].  For  the  sake 
of  completeness,  a  description  of  the  algorithm,  and  the  details  which  are  relevant 
to  its  hardware  implementation  are  given  below.  First,  we  review  key  definitions  as 
follows: 

Definition  1:  Let  Rp  be  an  image  row,  {\<p<m). 

(a)  A  run  in  Rp  is  a  sequence  of  1-pixels  (i.e.  pixels  with  gray  value  1)  of  Rp 
bounded  on  both  sides  by  0-pixel  (For  simplicity,  we  add  two  additional  0-pixels  to 
the  left  and  right-end  of  each  row  of  the  image,  respectively). 

(b)  The  lower  semi-image  Ip  consists  of  the  union  of  all  rows  Rj,  p^j^m. 

(c)  Gp  is  defined  to  be  the  partition  of  the  set  of  runs  in  Rp,  for  which  two  runs  are 
in  the  same  partition  group  g^Gp,  iff  they  belong  to  the  same  connected  com- 
ponent of  the  lower  semi-image  Ip. 

The  algorithm  consists  of  two  passes.  Pass  1  sweeps  through  the  rows  from  bot- 
tom to  top,  during  which  the  groups  in  Gp  are  calculated  inductively  from  the 
knowledge  of /?p  +  i  by  a  simple  updating  rule.  Pass  2  sweeps  from  top  to  bottom  to 
assign  component  numbers  to  each  pixel  and  outputs  this  symbolic  image. 

In  each  row,  runs  belonging  to  the  same  group  are  associated  by  a  simple 
mechanism,  called  bracket  marking,  as  follows: 

Definition  2:  We  define  four  symbols  ' [' ,  ']',  '[]'.  and  ']['. 

(a)  If  r  is  the  leftmost  (resp.  rightmost)  run  in  its  group  g,  it  is  assigned  marking  { 
(resp.  ]). 

(b)  If  r  is  both  leftmost  and  rightmost  run  in  g,  i.e.  g  consists  of  a  single  run  r, 
then  r  is  assigned  the  marking  [] . 

(c)  If  r  lies  between  the  leftmost  and  rightmost  run  in  its  group,  it  is  assigned  the 
marking  ][. 

An  example  is  given  in  Fig  2.1,  where  a  bracket  marking  is  shown  for  seven 
runs  in  Rp  numbered  from  right  to  left.  These  seven  runs  are  divided  into  three 
groups  -  (7,4,3,1),  (6,5),  and  (2).  It  is  shown  in  [8]  that  the  bracket  sequence  asso- 
ciated with  row  Rp  is  properly  nested,  i.e.  each  right  bracket  is  matched  to  an  asso- 
ciated left  bracket  and  vice  versa.  Hence,  the  groups  into  which  we  divided  the  set 
of  all  runs  in  Rp  can  be  reconstructed  from  the  bracket  sequence  associated  with  Rp. 

2.1  Pass  1 

The  goal  of  pass  1  is  to  calculate  a  bracket  marking  for  each  row  by  sweeping 
from  bottom  to  top,  and  to  store  it  in  a  bracket  table  (to  be  accessed  in  pass  2).  For 
this  we  want  a  rule  telling  us  how  to  calculate  the  bracket  marking  for  the  row  Rp 
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Figure  2.1  An  example  of  bracket  marking  assignment 

given  the  same  information  for  Rp  +  i.  The  key  observation  is  that  any  two  runs  in 
row  Rp  will  be  in  a  group  g^Gp  iff  there  are  some  runs  in  ^'€Gp+i  whose  index  set 
has  a  nonempty  intersection  with  the  index  sets  of  the  two  runs  in  Rp.  Intuitively, 
those  runs  in  row  Rp  which  are  linked  together  by  a  chain  in  the  semi-image  Ip 
belong  to  the  same  group.  The  algorithm  scans  two  rows,  Rp  and  Rp  +  i,  simultane- 
ously from  right  to  left.  Two  pushdown  stacks,  stackl ,  and  stack!,  are  used  in  pass 
1.  Stackl  is  used  for  tracing  the  group  relationship  represented  by  the  bracket 
sequence  in  Rp  +  i  as  following: 

(a)  If  a  run  with  ']'  is  encountered,  we  push  a  unmarked  symbol  onto  stackl; 

(b)  If  a  run  with  ']['  is  encountered,  we  neither  push  nor  pop; 

(c)  If  a  run  with  '['  is  encountered,  we  pop  the  top  entry  off  stackl  at  the  end  of  this 
run; 

(d)  If  a  run  with  '[]'  is  encountered,  we  push  an  unmarked  symbol  onto  stackl  at  the 
beginning  of  the  run,  and  pop  stackl  at  the  end  of  the  run. 

When  we  are  within  a  run  in  Rp  +  i  during  scan,  if  we  encounter  a  run  in  Rp,  then 
we  mark  the  symbol  on  the  top  of  stackl ,  which  indicates  that  the  current  group  in 
Rp  +  l  has  touched  some  run  in  Rp.  When  a  run  in  Rp  is  encountered  and  we  are 
within  a  run  in  Rp  +  i  with  a  ']['  or  '['  marking,  we  look  at  the  top  entry  of  stackl  to 
see  if  it  is  marked,  i.e.  to  see  if  the  current  group  in  Rp+i  already  touched  previ- 
ously some  run  in  Rp.  If  so,  the  current  encountered  run  in  Rp  is  linked  to  some  run 
to  its  right  in  the  same  row  through  a  path  below  Rp. 

Stackl  is  used  for  maintaining  the  group  information  about  the  row  Rp  during 
the  scan.  Assume  we  number  each  run  in  row  Rp  from  right  to  left  starting  from  1, 
as  an  example  shown  in  Fig.  2.1.  Each  entry  in  stackl  is  a  pair  of  integers, 
representing  the  right-most  and  left-most  runs  of  each  group  known  to  the  point  we 
have  reached  during  scan.  This  group  boundary  information  is  updated  according  to 
the  new  knowledge  we  gain  about  each  group  as  the  scan  continues.  Given  the  way 
that  we  scan  a  row  from  right  to  left,  the  connectivity  of  runs  in  row  Rp  through 
links  in  the  semi-image  /p  +  i  can  be  classified  into  three  basic  cases  as  follows  (refer 
to  Fig.  2.2): 

(a)  Several  runs  may  be  linked  by  a  simple  chain; 

(b)  A  cluster  of  links  (or  arcs)  may  be  bound  together  on  its  right  end  by  a  single 
run  in  row  Rp-, 

(c)  A  cluster  of  links  (or  arcs)  may  be  tied  together  on  its  left  end  by  a  single  run  in 
row  Rr,. 
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Figure  2.2  Basic  cases  of  connectivity  in  semi-image  Ip-. 

(a)  simple  chain;  (b)  a  cluster  of  links  is  bound  at  right  end; 

(c)  a  cluster  links  is  bound  at  left  end. 

The  general  case  can  be  decomposed  recursively  into  these  three  elementary 
cases  by  recursively  considering  a  subset  of  consecutive  runs  in  a  group  as  a  "single 
run".  For  example  in  Fig.  2.1,  if  we  consider  runs  4  and  3  as  a  "single  run"  denoted 
by  4-3,  the  connectivity  between  7,  4-3,  and  1  falls  into  the  elementary  case  (c). 
Our  algorithm  handles  the  connectivity  problem  from  inner  level  to  outer  level 
recursively.  Thus,  at  any  time,  we  only  have  to  deal  with  one  of  the  three  elemen- 
tary cases  defined  above;  however  a  run  could  be  a  simple  run  or  a  virtual  "single 
run". 

2.1.1  Manipulation  of  Stack  2  Information 

The  condition  for  beginning  a  new  group  can  be  easily  checked  locally.  When 
we  encounter  a  new  run  in  row  Rp,  we  check  if  it  makes  contact  with  any  run  in  row 
Rp  +  l.  If  it  does  not,  we  know  that  it  is  the  start  of  a  new  group  containing  one  run 
only.  On  the  other  hand,  if  it  does  make  contact  with  some  runs  in  Rp  +  i  the  right 
most  of  which  has  a  ']'  marking,  we  assume  temporarily  that  a  new  group  starts. 
(We  say  'temporarily'  here  since  initially  no  link  connects  the  current  run  in  Rp  to  a 
run  to  its  right  on  same  row.  However  this  new  run  might  link  to  some  run  to  its 
right  indirectly  through  several  arcs.  For  example  in  Fig.  2.2(c),  run  2  links 
indirectly  to  run  1  through  two  arcs.) 

Case  (a) 

In  case  (a)  of  Fig.  2.2,  we  have  entry  [1,1]  on  stack!  after  reaching  run  1,  which 
means  the  left  and  right  bound  of  the  current  group  are  1  and  1  respectively.  The 
left  bound  expands  to  2  after  run  2,  resulting  in  the  top  stack  entry  [1,1]  being 
updated  to  [2,1].  It  is  then  updated  to  [3,1]  when  run  3  is  reached,  after  which  it  is 
popped  off  stack!  since  a  group  is  finished.  So  in  case  (a),  the  group  always 
expands  its  left  bound. 
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Case  (b) 

The  situation  in  case  (b)  of  Fig.  2.2  is  slightly  more  complicated.  We  have  [1,1] 
on  stack!  after  run  1.  (Remember  that  we  have  two  marked  entries  on  stackl  after 
run  1  in  row  Rp  since  it  touched  two  runs  in  row  /?p  +  i  both  of  which  have  ']'  brack- 
ets.) 

After  run  2,  the  top  entry  on  stackl  becomes  [2,1]  and  the  top  entry  on  stackl 
was  popped.  At  this  point  we  do  not  know  whether  or  not  to  pop  the  top  entry  on 
stackl.  While  it  is  possible,  using  only  local  information  to  determine  the  start  of  a 
new  group,  this  is  not  the  case  for  determining  the  end  of  the  current  group.  The 
fact  that  a  run  has  no  more  arcs  linking  it  leftwards  does  not  necessarily  mean  that 
the  end  of  this  group  has  been  reached.  (Note  that  run  2  is  further  linked  to  run  3 
through  a  indirect  path.) 

To  handle  this  situation,  we  introduce  an  auxiliary  field  in  stackl 's  entry  to  indi- 
cate if  the  arc  associated  with  a  ']'  marking  in  that  entry  is  the  outermost  (or  lowest) 
arc.  For  example,  we  say  arc  A  is  the  lowest  among  the  cluster  of  arcs  bounded  by 
run  1  in  row  Rp  in  (b)  of  Fig. 2. 2.  The  auxiliary  bit  can  be  easily  set  as  following. 
During  a  run  in  row  Rp,  if  one  or  more  entries  have  to  be  pushed  onto  stackl,  the 
auxiliary  field  in  the  first  pushed  entry  is  set  to  TRUE  and  the  auxiliary  fields  in  all 
other  entries  are  set  to  FALSE.  If  the  auxiliary  field  in  the  entry  just  popped  off 
stackl  is  FALSE,  we  know  the  arc  (or  link)  just  ended  is  not  the  outermost  one 
(lowest),  i.e.  there  is  an  arc  belonging  to  the  same  group  enclosing  this  finished  arc, 
so  we  don't  pop  the  entry  off  stackl  because  we  don't  know  if  the  current  group  in 
row  Rp  is  finished  yet.  After  reaching  run  3  in  row  Rp,  the  left  bound  of  current 
group  is  expanded  to  3,  and  the  top  entry  on  stackl  is  changed  from  [2,1]  to  [3,1]. 
When  the  entry  on  stackl  is  popped  off,  we  see  that  the  auxiliary  field  is  TRUE.  So, 
we  can  confidently  pop  the  top  entry  on  stackl  off  this  time  because  we  know  that 
there  is  no  further  arc  linking  this  group  to  any  run  to  its  left. 

Case  (c) 

A  situation  symmetrical  to  case  (b)  occurs  in  case  (c).  After  reaching  run  1,  we 
have  one  entry,  [1,1],  on  stackl.  Note  that  at  the  moment  after  reaching  run  2,  it  is 
impossible  to  know  that  this  run  is  actually  linked  to  run  1  through  an  indirect  path. 
So,  we  have  a  second  entry,  [2,2],  on  top  of  stackl.  After  point  xi  in  run  3  is 
reached,  the  top  entry  of  stackl  is  updated  from  [2,2]  to  [3,2].  When  we  reach 
point  X2  we  find  that  two  arcs  are  bound  together.  So,  the  two  top  entries  on  stackl 
are  merged  into  one  in  such  a  way  that  the  new  left  bound  is  the  left  bound  of  first 
entry  on  stackl  and  the  new  right  bound  is  the  right  bound  of  the  second  entry  of 
stackl,  resulting  a  [3,1]  on  stackl.  These  operations  can  be  performed  in  two  steps: 
(1)  pop  stackl;  (2)  modify  the  left  bound  of  top  entry  by  the  left  bound  of  just 
popped  entry.  In  this  example,  we  pop  [3,2]  off  stackl  first;  then  change  the  left 
field  of  top  entry  -  [1,1]  to  3  and  get  [3,1]  as  the  new  top  on  stackl. 

2.1.2  Bracket  Table 


The  bracket  marking  can  be  effectively  encoded  by  a  two-bit  binary  number  as 
shown  in  Table  2.1.  The  bracket  information  calculated  during  pa^.?  1  is  stored  in  a 
memory  of  512x256x2  bits,  called  the  brackets  table.    This  table  is  indexed  by  two 
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Table  2.1  Bracket  Marking  Encode  Definitions 


Bracket  Marking 

Left  Bit 

Right  Bit 

[] 

0 

0 

[ 

0 

1 

] 

1 

0 

][ 

1 

1 

row  counters,  row_countl  and  row_count2,  and  two  column  counters,  run_countl  and 
run_count2 .  Row_countl  and  run_countl  are  combined  to  access  brackets  in  row  Rp, 
and  row_count2  and  run_count2  for  those  in  row  Rp  +  i.  Each  entry  in  the  table  has 
two  1-bit  fields,  both  of  which  are  initialized  to  zero.  Run_counts  are  set  to  zero  at 
the  beginning  of  a  row  sweep  and  incremented  by  1  every  time  a  new  run  started. 

Now  we  show  how  the  bracket  table  can  be  updated  efficiently  as  groups  grow 
during  scan.  There  are  only  two  types  of  group-expanding  operations:  (a)  the  left 
boundary  of  the  current  group  expands  to  include  the  current  scanned  run;  (b)  the 
two  top  most  groups  on  stack2  are  merged  together.  To  do  each  of  these  operations 
in  constant  time,  we  must  be  able  to  directly  index  correct  columns  in  a  row  of  the 
bracket  table.  The  index  of  the  current  scanning  run  is  provided  by  run_countl  for 
row  Rp  and  run_count2  for  row  Rp  +  \.  The  boundaries  of  all  non-ended  groups  are 
maintained  in  stack2  in  a  proper  nested  order  with  the  current  group  on  the  top. 
Thus,  the  updating  of  bracket  table's  entry  in  every  possible  case  can  be  done  in 
constant  time.  By  the  end  oi  pass  7,  we  get  a  complete  bracket  marking  for  every 
run  of  the  binary  image  stored  in  the  bracket  table. 

2.2  Pass  2 


In  pass  1  we  have  only  solved  the  connectivity  problem  for  runs  in  each  row  that 
connect  through  some  links  below  that  row.  However,  the  connectivity  between 
those  runs  in  a  row  which  are  linked  together  only  through  some  arcs  above  that 
row  can  not  be  found  in  this  bottom-to-top  sweep.  A  second  pass  sweeping  in  the 
opposite  direction  is  necessary  to  solve  the  connected  components  problem  com- 
pletely. 

Pass  2  is  much  simpler  than  pass  1 .  On  each  row,  two  sweeps  are  used.  The  first 
sweep  scans  the  pixels  in  two  consecutive  rows,  Rp-\  and  Rp,  simultaneously  from 
right  to  left.  The  second  sweep  scans  row  Rp  from  left  to  right  only.  We  still  use 
two  stacks,  stacks  and  stack4,  but  they  play  slightly  different  roles  to  the  two  stacks 
used  in  pass  1.  StackS  is  used  to  trace  the  group  information  in  row  Rp,  and  stack4 
saves  the  component  number  assigned  to  each  group  in  row  Rp.  A  circular  queue  is 
used  for  temporarily  storing  component  numbers  corresponding  to  each  run  in  row 
Rp-\  with  the  right-most  run  at  the  head  of  the  queue. 

The  first  sweep  induces  the  component  number  already  assigned  to  runs  in  row 
Rp-\  to  the  runs  in  row  Rp  to  which  they  are  connected,  and  allocates  a  new  com- 
ponent number  to  each  new  group  emerged  in  row  Rp. 

When  a  run  with  ']'  bracket  is  encountered  in  row  Rp,  which  means  the  start  of  a 
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new  group  has  been  found,  we  push  a  NULL  onto  stackS .  On  the  other  hand,  when 
a  run  is  encountered  in  row  Rp-\,  we  remove  the  entry  at  the  tail  of  the  queue 
which  contains  the  component  number  corresponding  to  this  run.  Then  if  the 
current  group  of  row  Rp  connects  with  any  run  in  row  Rp-\  and  the  top  of  stacks  is 
still  NULL,  we  replace  the  top  of  the  stacks  by  the  entry  most  recently  removed 
from  the  queue.  When  a  run  with  a  '['  bracket  is  encountered  in  row  Rp,  which 
indicates  the  current  group  has  ended,  we  pop  the  top  entry  from  stackS .  If  this 
popped  element  is  not  NULL,  it  is  pushed  onto  stack4;  otherwise,  a  new  component 
number  is  allocated  and  is  pushed  onto  stack4.  At  the  end  of  sweep  1,  we  have  the 
component  numbers  for  each  groups  of  row  Rp  in  stack4  with  the  last  ended  group's 
on  the  top. 

The  second  sweep  scans  row  Rp  only,  from  left  to  right.  When  a  run  with  '[' 
bracket  is  encountered,  we  pop  its  component  number  from  stack4  and  push  it  onto 
stacks .  The  component  number  of  a  group  on  stackS  is  popped  off  at  the  end  of  the 
group's  right-most  run  which  must  have  a  ']'  bracket.  Each  pixel  scanned  is  checked 
to  see  if  it  is  white.  If  it  is,  then  the  background  component  number  is  output.  Oth- 
erwise it  is  within  a  run,  so  we  output  the  content  on  the  top  of  stackS  which  is  the 
component  number  of  current  group. 

The  full  algorithm  written  in  C  Programming  Language  is  given  in  Appendix  A. 

3.  Hardware  Design 

In  section  2,  we  have  shown  that  all  operations  involved  in  the  computation  are 
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Figure  3.1  Block  Diagram  of  Connected  Component  Board 
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simple  and  can  be  classified  as  follows: 

(a)  Simple  variable  assignment  (register  transfer). 

(b)  Indexed  memory  access  (e.g.  bracket  table  read/write). 

(c)  Push  stack,  pop  stack,  or  modify  the  top  of  stack. 

(d)  Insert  queue  and  delete  queue. 

It  should,  and  does  indeed,  have  a  straightforward  implementation  in  hardware. 

The  goal  of  the  first  prototype  connected  components  board  designed  was  to 
implement  the  algorithm  described  in  section  2  in  standard  MSI  or  LSI  chips  without 
involving  special  VLSI  chip  design,  and  the  computation  on  a  512x512  binary  image 
can  be  accomplished  in  few  video-frame  times.  Real-time  implementation  and  VLSI 
design  are  discussed  in  the  last  section  of  this  paper. 

The  design  is  carried  out  in  a  systematic  way.  We  start  by  specifying  modular 
units  and  connections  between  them.  Then  the  data  structure  and  functions  per- 
formed in  each  module  are  precisely  described  by  a  number  of  processes  written  in 
a  special  microparallel  system  hardware  design  language  (see  Appendix  B). 
Processes  are  then  mapped  into  circuits  by  several  regular  rules.  Our  explanation 
concentrates  on  the  architecture  of  the  hardware.  The  low  level  schematic  design 
can  be  found  in  [10]. 

The  block  diagram  of  the  connected  component  board  is  shown  in  Fig.  3.1. 
There  are  five  major  components: 

(1)  Pass  1  processor. 

(2)  Pass  2  processor. 

(3)  Image  buffer. 

(4)  Bracket  table. 

(5)  Input/Output  interface. 

In  following,  we  examine  each  unit  by  describing  the  main  data  structure  of  the 
component  along  with  its  control  structure. 

3.1  Control  Units 


The  computation  of  the  algorithm  is  carried  out  by  Pass  1  Processor  and  Pass  2 
Processor.  Both  processors  have  very  similar  both  data  and  control  structures. 

As  can  been  seen  in  Fig.  3.1,  Pass  1  processor  contains  a  control  unit,  two 
stacks,  and  some  local  variables.  The  control  unit  coordinates  the  register  transfers, 
stack  operations,  image  buffer  and  bracket  table  accessings,  according  to  the  algo- 
rithm by  issuing  pulses  to  each  component  in  an  appropriate  order.  The  control  unit 
is  implemented  by  a  typical  time  sequential  circuit  or  so  called  finite  state  machine. 
The  structure  of  a  finite  state  machine  is  given  in  Fig.  3.2  and  its  principle  is 
explained  in  many  digital  logic  design  books,  e.g.  [11].  The  output  signals  are  con- 
trol signals  to  each  data  components.  The  logical  structure  of  control  flow  in  the 
control  unit  of  pass  1  is  depicted  in  Fig. 3. 3. 

There  are  two  nested  loops  in  the  diagram,  with  the  outer  one  iterating  on  rows 
and  the  inner  one  on  columns.  There  are  number  of  modules,  so  called  processes,  in 
Fig.  3.3.  All  processes  are  self-timed.  Within  each  process  some  control  signals  are 
issued  to  data  components  to  carry  out  manipulations.  For  example,  when 
Passljnitialization  process  is  activated,  it  sends  control  signals  to  set  initial  state  for 
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Figure  3.2  Block  diagram  of  finite  state  machine 
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Figure  3.3  The  control  unit  of  Pass  1  processor 

pass  1  in  various  data  components.    The  procedural  descriptions  of  those  processes 
can  be  found  in  Appendix  B. 

The  control  unit  of  the  Pass  2  Processor  has  a  similar  structure  as  the  Pass  1 
Processor.  Fig.  3.4  shows  the  logical  control  structure  of  the  Pass  2  control  unit.  It 
contains  three  loops.  The  outer  loop  iterates  on  rows.  The  two  inner  loops  are  con- 
catenated in  serial.  The  first  inner  loop  iterates  on  columns  from  right  to  left  in  a 
row;  the  second  loop  iterates  in  the  opposite  direction. 

3.2  Data  Components 

3.2.1  Pushdown  Stacks 

The  most  important  data  structure  here  is  the  pushdown  stack.  The  depth  of 
each  stack  used  is  256.  A  bi-directional  shift-register  is  seems  the  best  element  for 
implementing  a  pushdown  stack.  Since  we  couldn't  find  commercial  shift-registers  of 
such   length,   each   stack   in   our   design   is   constructed   painfully   from  a  piece  of 
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memory  and  a  counter  serving  as  the  stack  pointer  as  depicted  in  Fig.  3.5.  Since  the 
content  of  the  top  of  a  stack  is  frequently  referenced  as  variables  throughout  the 
process,  it  is  stored  in  a  register  rather  than  in  memory.  A  push  stack  operation  is 
performed  in  three  steps: 
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(1)  Increment  stack  pointer  by  one. 

(2)  Write  the  content  of  the  register  (the  last  top  entry)  into  memory. 

(3)  Gate  the  external  value  (the  new  top  entry)  into  the  register. 

The  pop  operation  is  done  in  a  similar  manner.  The  updating  of  the  top  entry  of 
a  stack  can  be  most  efficiently  done  by  just  gating  the  external  input  to  the  register. 

Note  that  since  the  stack  operation  needs  more  than  one  clock  cycle,  we  have  to 
latch  the  external  input  and  control  signals  to  the  stack  during  the  whole  stack 
operation.  Such  latching  circuits  can  be  avoided  if  we  design  a  stack  in  VLSI  by  just 
a  simple  array  of  bi-directional  shift-register,  since  any  operation  needs  only  one 
cycle  in  this  case. 

3.2.2  Image  Buffer 

The  image  buffer  is  for  storing  a  512x512  1-bit  binary  image  whose  structure  is 
shown  in  Fig.  3.6.  Because  of  the  relative  low  speed  of  memory  access,  we  compact 
8  bits  into  one  word.  So  we  need  a  total  of  32K  bytes  memory.  While  reading  pix- 
els, a  word  is  loaded  from  the  memory  into  an  8-bits  shift-register.  Then  individual 
bits  can  be  get  at  the  output  pin  of  the  shift-register  sequentially.  The  advantage  of 
this  is  that  we  have  only  one  memory  access  for  every  eight  pixels.  In  addition,  the 
loading  of  a  word  from  memory  to  shift-register  could  be  overlapped  with  pixel 
computation. 

As  we  mentioned  in  the  algorithm,  we  need  to  scan  two  rows  simultaneously. 
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Since  the  two  rows  that  we  want  to  read  at  the  same  time  are  always  consecutive, 
one  is  an  even  numbered  row  and  the  other  is  odd  numbered.  This  suggests  a 
natural  division  of  the  image  buffer  into  two  ranks,  one  for  even  rows  and  one  for 
odd  rows,  so  we  don't  have  conflicts  in  reading  two  consecutive  rows  at  a  time.  The 
pixels  from  the  current  row  and  the  previous  row  are  obtained  at  out  pins  of  shift- 
registers  and  latched  in  D-ff's. 

Row_counterl,  row_counter2,  and  column_counter  provide  the  effective 
addresses  for  accessing  two  memory  ranks  through  an  address  decoder. 

3.2.3  Bracket  Table 

The  bracket  table  is  shown  in  Fig.  3.7.  The  bracket  markings  are  encoded  by 
two  bits.  These  two  bits  are  stored  in  two  ranks  of  memory,  512x256x1  bits  each, 
for  left-  and  right-bits  respectively.  The  current  bracket  markings,  read  from  the 
table,  are  latched  in  two  D-ffs. 

The  row  number  is  provided  by  row_counterl  and  row_counter2  through  a  mul- 
tiplexer. The  column  number  is  chosen  from  run_counterl,  run_counter2,  the  top  of 
stack!  or  last  popped  entry  from  stack! .  An  address  decoder  generates  effective 
addresses  for  both  ranks. 

3.3  Input/Output  Interface 
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The  general  purpose  handshaking  interface  we  designed  is  flexible  enough  to  be 
attached  to  a  DMA  controller,  a  parallel  port  or  a  serial  port  easily  with  a  little 
additional  logic. 

The  block  diagram  of  the  interface  is  shown  in  Fig.  3.8.  There  are  two  registers: 
(1)  a  status  register;  and  (2)  a  command  register.  The  meaning  of  each  bit  in  the 
status  register  is  given  in  Fig.  3.9(a)  and  in  the  command  register  in  Fig.  3.9(b). 
The  data  bus  has  eight  bit  bi-directional  lines.  In  the  write-to-board  mode,  the  data 
is  supplied  by  an  external  device;  and  in  read-from-board  mode,  the  data  bus  is  set 
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Figure  3.8  The  block  diagram  of  interface 
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by  the  board.  The  data  transfer  direction  is  controlled  by  two  input  lines,  PREAD 
and  PWRITE,  which  are  low  active  and  exclusive.  Another  output  line  from  the 
board,  called  READY,  is  an  acknowledge  signal  to  an  external  device  to  indicate 
that  the  buffer  register  of  the  board  is  empty  in  write-to-board  mode,  and  full  in 
read-from-board  mode.  The  signal  READY  is  high  active. 

The  timing  of  those  signals  for  I/O  operations  is  given  in  [10]. 

4.  Considerations  on  VLSI  design 
and  Real-time  Applications 

Since  the  connected  components  algorithm  is  such  a  fundamental  procedure  in 
image  analysis  and  has  more  and  more  applications  in  robot  vision,  it  is  very  attrac- 
tive to  implement  this  design  on  a  VLSI  chip,  so  it  can  be  conveniently  added  to  an 
image  processing  computer  or  general  purpose  system.  Although  we  can  simply 
map  the  whole  circuit  onto  a  chip  and  replace  each  logical  gate  by  a  corresponding 
cell  of  VLSI  technique,  it  is  not  the  best  way  to  do  it.  The  disadvantage  of  this 
approach  is  that  the  irregular  structure  of  original  design  makes  the  physical  layout 
in  VLSI  difficult.  Many  of  the  components  can  be  implemented  instead  in  a  more 
efficient  and  much  simpler  style  of  VLSI. 

For  instance,  a  pushdown  stack  can  be  implemented  by  an  array  of  bi-directional 
shift  registers.  It  has  several  advantages:  (a)  it  has  a  very  regular  pattern,  so  it  is 
extremely  simple  to  change  the  dimension  of  a  stack;  (b)  a  stack  pointer  is  no  longer 
necessary;  (c)  each  of  the  three  kinds  of  operations  on  a  stack  can  be  performed  in 
just  one  cycle;  (d)  because  of  (c),  the  latching  circuit  for  external  inputs  and  control 
signal  to  the  stack  can  be  eliminated. 

A  control  unit  in  VLSI  is  usually  implemented  in  a  more  regular  structure  rather 
than  a  random  logic  design.    The  typical  structure  of  a  finite  state  machine  in  VLSI 
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is  shown  in  Fig.  4.1.  The  combinational  logic  network  implements  logical  functions 
in  the  form  of 

;  =  i 
It  is  not  very  difficult  to  write  down  the  input  to  a  register  in  our  design  as  the  func- 
tion of  other  variables  in  the  above  form.  So  the  full  VLSI  implementation  is  obvi- 
ously possible. 

More  parallelism  can  be  obtained  at  both  the  global  and  local  level.  For  exam- 
ple, in  the  above  design  as  Pass  1  processor  is  running.  Pass  2  processor  is  idle. 
Globally,  it  is  valuable  to  pipeline  the  two  processors  when  we  need  to  compute  the 
input  images  continuously.  A  possible  configuration  is  shown  in  Fig.  4.2.  Instead  of 
putting  whole  system  on  one  chip,  it  can  be  naturally  divided  into  four  modules, 
namely  the  pass  1  processor,  the  pass  2  processor,  the  image  buffer,  and  the  bracket 
table.  In  Fig.  4.2  we  have  two  image  buffers  and  two  bracket  tables.  When  the  pass 
1  processor  is  working  on  set  1,  the  pass  2  processor  works  on  set  2.  At  the  end  of  a 
pass,  the  two  processors  can  exchange  image  buffers  and  bracket  tables.  So  when 
the  pass  1  processor  finishes  one  pass  on  image  1,  the  result  is  stored  in  bracket 
table  1.  Then  these  data  are  passed  to  the  pass  2  processor  and  at  the  same  time  the 
pass  1  processor  starts  working  on  the  next  incoming  image  with  the  other  image 
buffer  and  bracket  table.   Data  switching  can  be  done  trivially  with  multiplexers. 

If  the  order  of  output  sequence  of  a  symbolic  image  is  allowed  to  be  reversed, 
we  can  reverse  the  sweeping  directions  of  both  pass  1  and  pass  2.  In  such  a  case, 
the  pass  1  processor  doesn't  have  to  wait  until  the  entire  image  has  been  loaded  into 
the  image  buffer.  In  other  words,  it  can  run  during  the  image  loading  process.  This, 
however,  introduces  the  possibihty  of  simultaneous  accesses  to  the  image  buffer  by 
the  image  loading  process  and  the  pass  1  computation  process.  So  some  additional 
mechanism,  such  as  a  semaphore,  is  needed  to  coordinate  the  processes. 

Locally,  some  operations  in  each  of  the  processes  can  be  obviously  performed 
during  a  single  cycle  without  conflicts  in  the  data  path.    Such  parallelism  can  be 
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CCB  page  16 

easily  achieved  by  broadcasting  a  control  signal  to  more  than  one  proper  data  com- 
ponents. This  kind  of  parallelism  has  been  realized  in  the  many  processes  of  our 
design. 

As  the  speed  of  computation  increases,  data  transfer  will  cause  a  bottleneck.  It 
seems  that  regular  DMA  transfer  in  micro  or  mini  systems  would  not  be  able  to 
match  the  the  computation  speed.  An  ultra  fast  secondary  bus  directly  connecting  it 
to  other  system  modulars  is  more  preferable. 
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Connected  Components  Algorithm 
Implementated  in  C  Programming  Language 


/"  Constants  and  macro  definitions  */ 


#lnclude 
#include 


<stdio.h> 
<math.h> 


#deflne  MAXROW 
#deflne  MAXCOL 
#deflne  MAXRUN 
#deflne  Qlcn 

#define  START_C 
#deflne  START_P 
#deflne  END_C 
#define  END_P 
#deflne  OVERLAP 
#denne  IN_PRUN 
#deflne  IN_CRUN 

#deflne  BRKT 
#denne  LBRKT 
#define  LEFT 
#define  RIGHT 
#deflne  CLOSE 
#deflne  MID 


512 

512 

(MAXCOL/2) 

(2*MAXRUN) 

(!(l_cpix)&&(cpix)) 

(!(Lppix)&&(ppix)) 

(!(cput)&&(l_cpix)) 

(!(ppix)&&(l_ppix)) 

((cpu)&&(ppix)) 

(ppix) 

(cpix) 

brkt[01*2  +  brkt[l] 

Lbrkt[0]*2  +  Lbrkt[l] 

1 

2 

0 


/'  Cloabt  subroutine  */ 
Int  conncomps(); 

/*  Global  Variables   '/ 

short       BrktL[MAXROW][MAXRUN];  I*  left  columun  of  brackets  table' I 

short       BrktR[MAXROW][MAXRUN];  /*  right  columun  of  breackets  table  'I 

Int  run_COuntl;  /•  current  run's  number  In  current  row  */ 

int  run_count2;  /*  current  run's  number  in  previous  row  */ 

Int  brkt[2];  /•  bracket  marking  associated  with 

the  current  run  in  previous  row  */ 

int  Lbrkt[2];  I*  last  content  of  brkt[)  '! 


extern     short       image[MAXROW  +  2][MAXCOL+2];  I*  source  image  'I 

extern     short       components[MAXROW][MAXCOL];  /*  destination  image  ' 

Int  colour[7]  /•  codes  for  VICOM  pseudo  colours  'I 

=  {0X120,  0X5A6,  0X240,  0X364,  0X6C8,  0X000,  0X486}; 
int  cpix,  /*  current  pixel  In  current  row  */ 

ppix,  /*  current  pixel  in  previous  row  */ 

Lcpix,  /•  last  pixel  in  current  row  'I 

l-ppix;  /*  last  pixel  in  previous  row  'I 
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int  queue_countl,queue_count2; 

Int  componcnt_No; 


/'  components  number  queue  pointers  */ 
/*  next  new  identify  number  */ 


/*  Main  subroutine  */ 


connected_conncomps() 
{ 


int 


row, col; 


/•  initial  brackets  table  */ 

for  (row  =  0;  row  <  MAXROW;  row+ +) 

for  (col  =  0;col  <  MAXRUN;  col+ +)  { 

BrktL[row][col]  =  BrktR[row][col]  =  0; 

} 

/*  PASS  1:  Scan  from  bottom  to  top  to  generate 

barckets  marking  for  each  run  of  each  row  */ 

for  (row  =  MAXROW;row>0;row--) 
Swcep_of_Passl(row) ; 

/*  PASS  2:  Sacn  from  top  to  bottom  to  assign 

components  ID#  to  each  pixel  of  the  Image  */ 

component_No  =  1; 
queue_countl  =  queue_count2  =  0; 
for  (row  =  l;  row<  =  MAXROW;  row+  +  ) 
Swcep_of_Pass2(row); 
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Scan  row  for  PASS  J  to  ger\erale  brackets  marking 
for  each  run  of  the  row  scanned. 


•/ 


Sweep_of_Passl(i) 

int  i; 

{ 


short       stackl[MAXRUN][2]; 


int 
Int 


spl; 
poppedl; 


short       stack2[MAXRUN][3]; 


Int 

sp3; 

int 

poppcd2[3]; 

int 

touch_p; 

int 

touch_c; 

int 

col; 

/*  current  scanned  row  number  */ 

/*  stackl  Is  for  tracking  brackets  in  previous  row. 
stackl[][0]  =  1 ,  Indicates  the  current  group  in 
previous  row  contacted  one  or  more  runs 
in  current  row, 
=  0,  otherwise. 

stackl  [III],  auxiliary  flag.  */ 
/*  stack  pointer  for  stackl  */ 
/•  last  content  of  stackl[spl] [1]  'I 

/*  stackl  Is  for  malntanlng  the  temporary  boundaries 

of  groups  in  current  row. 

stack2[][0],  right-most  run  of  current  group; 

Stack2[] [1] ,  left-most  run  of  current  group; 

stack2f]f2],  auxiliary  flag.  */ 
/*  stack  pointer  for  stack2  */ 
/*  last  entry  popped  from  stack2  */ 

/*  =1,  Indicates  the  current  run  in  current 

row  already  touched  an  run  in  previous  rowl  'I 
/*  current  run  In  previous  row  already  touched 
an  run  In  current  row  */ 


/*  Initialization  at  beginning  of  a  sweep  */ 

touch_p  =  touch_c  =  0; 
spl  =    sp3  =  0; 
cpix  =  ppix  =  0; 
run_countl  =  run_count2  =  0; 

for  (col=MAXCOL;col  >=0;  col--)  { 
Lcpix  =  cpix; 
Lppix  =  ppijt; 
cpix  =  imagc[i][col]; 
ppix  =  image[i+l][coI]; 


/*  backup  last  pixel  *l 

/•  current  pixel  In  current  row  '/ 
/'  current  pixel  In  previous  row  */ 


/*  if  an  run  starts  in  previous  row,  read  the  associated  bracket 

marking  from  the  table    */ 
if  (START_P)  { 

Lbrkt[0]  =  brkt[0]; 

Lbrkt[l]  =  brkt[l]; 

run_count2+  + ; 

brkt[0]  =  BrktL[i+l][run_count2]; 

brkt[l]  =  BrktR[i+l][run_count2]; 
} 
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/*  if  an  run  starts  In  current  row,  initial  the  entry  of  it  in  the  table  */ 
If  (START.C)  { 

run_countl+  +  ; 

BrktL[i][run_countl]  =  0; 

BrktR[i][run_countl]  =  0; 

} 

/*  runs  start  in  both  current  row  and  previous  raw  simultaneously  */ 
If  ((START_P)  &&  (START.C))  { 
switch  (BRKT)  { 
/*  push  a  marked  entry  on  stack  1 ; 

push  a  new  group  entry  on  stack2.  */ 
case  RIGHT: 
case  CLOSE: 

spl+  +  ; 

stackl[spl][0]  =  stackl[spl][l]  =  1; 

sp3+  +  ; 

stack2[sp3][0]  =  stack2[sp3][l]  =  run_countl; 

break; 

case  MID    : 
case  LEFT  : 

/*  expand  the  group  on  the  top  of  stack!  at  its  right  boundary  */ 
If  (suckl[spl][0])  { 

BrktL[i][siack2[sp3][0]]  =  1; 
BrktR[i][run_countl]  =  1; 
stack2[sp3][0]  =  ruii_countl; 

} 

/*  push  a  new  group  entry  on  slack!  */ 
else  { 

sp3+  +  ; 

stack2[sp31[0]  =  stack2[sp3][l]  =  run.countl; 

If  (BRKT  =  =  MID)  { 

$tackl[spl][0]  =  1; 

} 
} 
break; 

} 
} 

/*  an  run  starts  in  previous  row,  but  no  transition  in  current  row  'I 
else  If  ((START_P)  &&  !(START_C)  &&  !(END_C))  { 
If  (lOVERLAP)  { 

/*  push  an  unmarked  entry  on  stack!  */ 
•witch  (BRKT)  { 
case  RIGHT: 
case  CLOSE: 

spl+  +  ; 

stackl[spl][0]  =  0; 
stackl[spl][l]  =  1; 
break; 

} 
}  else  { 
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switch  (BRKT)  { 

/*  push  a  marked  entry  on  stackl  */ 
case  RIGHT: 
case  CLOSE: 

spH-  +  ; 

stackl[spl][0]  =  1; 
If  (!touch_p) 

$tackl[spl][l]  =  1; 
else  { 

If  (l_brkt[0]  =  =  1) 

stackl[spl][l]  =  0; 
else 

stackl[spl][l]  =  stack2[sp3J[2]; 
}  break; 
/*  Combine  two  top  groups  of  stackl  into  one  */ 
case  MID    : 

case  LEFT  :  If  (stackl[spl][0]  &&  ((touch.p  && 
!l_brkt[0]  &&  poppcdl)  || 
!touch_p))  { 

popped2[0]  =  stack2[sp3][0]; 
poppcd2[l]  =  stack2[sp3][l]; 
sp3--; 

BrktR[i][popped2[l]]  =  1; 
BrktL[i][stack2[sp3][0]]  =  1; 
stack2[sp3][0]  =  popped2[0]; 
}  break; 
} 
} 
} 

/*  an  run  starts  in  current  row,  but  no  transition  in  previous  row  */ 
else  If  ((START.C)  &&  !(START_P)  &&  !(END_P))  { 
If  (lOVERLAP)  { 

/*  push  a  new  group  entry  on  stack2  */ 

sp3+  +  ; 

stack2[sp3][0]  =  stack2[sp3][l]  =  run_countl; 

stack2[sp3][2]  =  1; 

}  el«  { 

switch  (BRKT)  { 
case  RIGHT: 
case  CLOSE: 

/*  mark  the  entry  on  the  top  of  stackl ; 

push  a  new  group  entry  on  stack2 .  */ 
If  (!touch_c)  { 

sp3+  +  ; 

stack2[sp3][0]  = 

stack2[sp3][l]  =  run_countl; 

stack2[sp3I[2]  =  1; 

stackl[spl][0]  =  1; 

stackl[spl][l]  =  1; 

} 

/*  expand  the  group  on  the  top  of  stackl  at  its  right  boundary  'I 

else  { 
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BrktR[i][run_countl]  =  1; 

BrktL[i][stack2[sp3][0]l=l; 
stack2[sp3][0]  =  run_countl; 

} 

break; 
case  MID    : 
case  LEFT  : 

/*  expand  the  group  on  the  top  of  stack!  at  its  right  boundary  */ 
If  (touch_c)  { 

BrktR[il[run_countl]  =  1; 
BrktL[i][stack2[sp3][0]]  =  l; 
stack2[sp3][0]  =  run_countl; 

}  «1»«  { 

/*  Expand  the  group  on  the  top  of  stack2  at 

lis  right  bundary  '/ 
If  (stackl[spl][0])  { 

BrktR[i][run_countl]=  1; 
BrktL[i][stack2[sp3][0]]  =  l; 
stack2[sp3][0]  =  run_countl; 

} 

/*  push  a  new  group  entry  on  stack!  */ 

cUc{ 

sp3+  +  ; 

stack2[sp3][0]  =  run_countl; 

stack2[sp3][l]  =  run_countl; 

If  (BRKT  ==  MID){ 

stackl[spl][0]=l; 

} 
} 
} 


} 


/•  both  runs  in  current  row  and  previous  row  terminate  simultaneously  */ 
else  If  ((END_P)  &&  (END_C))  { 
touch_c  =  touch_p  =  0; 

/•  If  the  group  on  the  top  of  stack!  is  finished,  pop  It  off  stack!.  'I 
switch  (BRKT)  { 
case  CLOSE: 
case  LEFT  : 

If  (stackl[spl][l])  sp3--; 

poppcdl  =  stackl[spl--][l]; 

break; 
} 
} 

/•  the  run  In  previous  row  terminates,  but  no  transition  In  current  row.  '! 
else  If  ((END_P)  &&  I(START_C)  &&  !(END_C))  { 
If  (!IN_CRUN)  { 

/•  pop  the  group  on  the  top  of  stack!,  If  it  Is  end.  '! 

switch  (BRKT)  { 
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case  CLOSE: 

if  ((touch_c)  &&  (stackl[spl][l])) 

sp3--; 
poppedl  =  stackl[spl--][l]; 
break; 
case  LEFT  : 

if  ((stackl[spl][0]  &&  stackl[spl][l]) 
II  (touch_c  &&  stack2[sp3][2])) 
sp3--; 
poppedl  =  stackl[spl--][l]; 
break; 
} 
}  el»e  { 

touch_p  =  1; 

switch  (BRKT)  { 

/*  backup  the  auxiliary  flag,  then  pop  stack] .  */ 

case  CLOSE: 

case  LEFT  : 

stack2[sp3][2]  =  stackl[spl][l]; 
poppedl  =  stackl[spl--][l); 
break; 

} 
}  touch_c  =  0; 

} 

/*  the  run  in  current  row  terminates,  but  no  transition  In  previous  row.  '/ 
else  if  ((END_C)  &&  !(START_P)  &&  !(END_P))  { 

/*  If  the  group  on  the  top  of  stack2  is  end,  pop  it  off  stackZ .  '! 
if  (!IN_PRUN)  { 

if  (!touch_p)  { 
sp3--; 
}  «lse  { 

switch  (BRKT)  { 
case  CLOSE: 

case  LEFT  :  If  (stack2[sp3][2])  { 

sp3--; 
}  break; 
} 
} 
}  else  { 

touch_c  =  1; 
}  touch_p  =  0; 
} 

/*  an  run  starts  In  previous  row,  and  the  run  in  current  row  terminates  at  the  same  time.  */ 
else  if  ((START_P)  &&  (END_C))  { 

/*  if  the  group  on  the  top  of  stack2  is  end,  pop  it  off.  'I 
if  (!touch_p)  { 
sp3--; 
}  else  ( 

switch  (BRKT)  { 

case  CLOSE: 

case  LEFT  :  if  (stack2[sp3][2])  { 
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sp3--; 
}  break; 

} 
} 

iwltch  (BRKT)  { 

/*  push  an  unmarked  entry  on  stack!  *l 
case  RIGHT: 
case  CLOSE: 

spl+  +  ; 

stackl[spl][0]  =  0; 

$tackl[spl][l]  =  1; 

break; 

} 

touch_p  =  touch_c  =  0; 

} 

/*  an  run  starts  in  current  row,  and  the  run  in  previous  row  terminates  at  the  same  time.  */ 
else  if  ((START.C)  &&  (END_P))  { 

switch  (BRKT)  { 

/*  if  the  group  on  the  top  of  stack2  Is  end,  pop  it  off.  */ 

case  CLOSE: 

case  LEFT  :         if  ((touch_c)  &&  (stackl[spl][l]))  { 
sp3-; 

} 

poppedl  =  stackl[spl--][l]; 

break; 

} 

/*  push  a  new  group  entry  o  slack2  */ 
sp3+  +  ; 

stack2[sp3][0]  =  $tack2[sp3][l]  =  run_countl; 
stack2[sp3][21  =  1; 
touch_p  =  touch_c  =  0; 
} 
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Swc€p_of_Pass2(i) 
Int  i; 


/•  stack!  is  for  tracking  brackets  in  previous  row  */ 
/'  slack  pointer  for  stack!  '/ 

/*  stack2  is  for  storing  components  numbers 
of  runs  in  current  row.  */ 


short       stackl[MAXRUN][2]; 
Int  spl; 

int  stack2[MAXRUN]; 

int  sp2; 

int  coniponcnts_qucue[QIcn]; 

int  col, code; 

Int  current_component; 

spl  =  sp2  =  0; 

cpix  =  ppix  =  0; 

run_countl  =  0; 

for  (col=MAXCOL;col  >=0;col--)  { 

Lcpix  =  cpix; 

Lppix  =  ppix; 

cpix  =  image[i][col]; 

ppix  =  image[i-l][col]; 

/*  if  an  run  strarts  in  previous  row,  increment  queue_count2 

to  point  to  the  entry  of  queue  containning  the  component  number 
associated  with  this  run.  */ 

If  (START_P)  { 

queue_count2  =  ((qucue_count2+ 1)  %  (Qlen)); 

/•  if  an  run  starts  in  current  row,  increment  queue_countl  to 

allocate  an  entry  in  queue  for  this  run.  Read  the  bracket  marking 
associated  with  this  run  from  the  table.  'I 
If  (START.C)  { 

queue_countl  =  ((queue_countl+ 1)  %  (Qlcn)); 
run_countl+  + ; 

brl£t[0]  =  BrktL[i][run_countl]; 
brkt[l]  =  BrktR[i][run_countl]; 
} 

If  ((START_P)  &&  (START.C))  { 

/•  induce  the  component  number  from  the  contacted  run  in  previous  row  */ 
iwitch  (BRKT)  { 
case  RIGHT: 
case  CLOSE: 

spl+  +  ; 

stackl[spl][0]  =  components  _queuc[queuc_count2]; 

$tackl[spl][l]  =  1; 

break; 
case  LEFT  : 
case  MID   :  if  (!stackl[spl][l])  { 

stackl[spl][0]  =  coniponents_queue[qucue_count2]; 
stackl[spl][l]  =  1; 
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} 

break; 

} 
}  else  if  ((START_P)  &&  !(START_C)  &&  !(END_C))  { 

/*  Induce  the  component  number  from  the  contacted  run  in  previous  row  */ 
if  ((OVERLAP)  &&  (!stackl[spl][l]))  { 

stackl[spl][0]  =  components_queue[queue_count2]; 

stackl[spl][l]  =  1; 

} 
}  else  if  ((START.C)  &&  !(START_P)  &&  !(END_P))  { 
switch  (BRKT)  { 
case  RIGHT: 
case  CLOSE: 

/*  set  the  component  number  to  NULL  temporarily  */ 
If  (lOVERLAP)  { 
spl+  +  ; 
stackl[spl][0]  =  stackl[spl][l]  =  0; 

} 

/*  induce  the  component  number  from  previous  row  */ 

else  { 

spl+  +  ; 

stackl[spl][0]  =  componcnts_queuc[qucuc_coun:2]; 

stackl[spl][l]  =  1; 
}  break; 
case  MID    : 
case  LEFT  : 

/*  induce  the  component  number  from  previous  row  */ 
If  ((OVERLAP)  &&  (!stackl[spl][l]))  { 

stackl[spl][0]  =  components_queuc[qucue_count2]; 

stackl[spl][ll  =   1; 
}  break; 

} 
}  else  if  ((END_P)  &&  (END.C))  { 

/*  save  the  component  number  for  a  finished  group  in  current 

row  into  stack!  to  be  used  in  second  sweep.  */ 
•witch  (BRKT)  { 
case  CLOSE: 
case  LEFT  : 

stack2[++sp2]  =  stackl[spl--][0]; 
break; 
} 
}  else  if  ((END_P)  &&  !(START_C)  &&  !(END_C))  { 

/•  NULL  •/ 
}  else  if  (((END.C)  &&  !(START_P)  &&  !(END_P))  || 
((START.P)  &&  (END_C)))  { 
switch  (BRKT)  { 
case  CLOSE: 
case  LEFT  : 

/•  save  the  component  number  of  a  finished  group  into  stack!  */ 
if  (stackl[spl][l]) 

stack2[++sp2]  =  stackl[spl--][0]; 
/*  if  the  finished  group  belongs  to  a  new  component,  then 
asslgne  a  new  number  to  It  and  save  in  stack!.  */ 
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elsc{ 

stack2[+ +sp2]  =  coinponent_No+ +  ; 
spl--; 
}  break; 

} 
}  else  If  ((START.C)  &&  (END_P))  { 

/*  set  the  component  number  to  NULL  temporarily  */ 
switch  (BRKT)  { 
case  RIGHT: 
case  CLOSE: 

spl+  +  ; 

stackl[spl][0]  =  stackl[spl][l]  =  0; 
break; 
} 
} 
} 

spl  =  0;cpix  =  0; 
currcnt_component  =  0; 
queue_count2  =  queue_countl; 

for  (col  =  0;col<  =  MAXCOL;coH-+)  { 
l_cpix  =  cpix; 
cpix  =  image[i][col]; 

If  (START_C)  { 

brkt[0]  =  BrktL[i][run_countl]; 
brkt[l]  =  BrktR[i)[run_countl]; 
run_countl--; 

switch  (BRKT)  { 

case  LEFT  :  current_componcnt  =  stack2[sp2--]; 

stackl[+ +spl)[0]  =  current_component; 

break; 
case  MID    :  current_componcnt  =  $tackl[spl][0]; 

break; 
case  RIGHT:        current_component  =  stackl[spl--][0]; 

break; 
case  CLOSE:        currcnt_component  =  stack2[sp2--]; 

break; 

} 

/'  save  the  component  number  of  current  run  Into  queue  */ 
components_queue[qucue_count2]  =  current_component; 
queue_coun:2  =  ((queuc_count2-l  +  Qlcn)  %  (Qlen)); 
}  else  if  (END_C)  { 

currcnt_componcnt  =  0; 

} 

/*  map  the  component  number  into  pseudo  color,  and  store  it 

into  the  destination  image.  */ 
code  =  (current_componcnt  =  =  0)  ?  6  :  (c_id  %  6); 
componcnts[i][col)  =  colour[code]; 
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Appendix.  B 


Process  Description  of  The  Algorithm  in 
A  Microparallel  Hardware  Design  Language 


B.l  Definitions  of  Variables  and  Flags 


image_buffer(0.. 511,0.. 511); 


Left_bracket_table(0.. 511,0.. 255); 
Right_bracket_table(0.. 511,0.. 255); 

Row_Countl; 

Row_Count2; 

Column_Count; 

Bracket_Countl; 

Bracket_Count2; 


Pl_Stackl(256,2) 

Pl_Stack2(256,3) 

P2_Stackl(256,2) 

P2_Stack2(256,l) 

Sl_popped(2); 

S2_popped(3); 

Queue(0..512); 
Queue_Countl; 
Queue_Count2; 

CPIX; 
LCPIX; 
PPIX; 
LPPIX; 

brkt(2); 
l_brkt(2); 

touch_c; 
touch_p; 

In_Down_load; 

In_Passl; 

In_Pass2; 


Start_passl; 
Start_pass2; 
Sweep2_of_Pass2_enable; 


/•  The  image  storage  buffer  is  divided 
into  two  ranks,  one  for  even 
rows  and  one  for  odd  */ 

*  Table  of  left  bracket  bits  */ 

*  Table  of  right  bracket  bits  */ 

*  9  bits:  binary  up-down  counter  *  I 

*  9  bits:  Binary  up-down  counter  *l 

*  9  bits:  binary  up-down  counter  *l 

*  8  bits:  up -down  counter  */ 

*  8  bits:  up -down  counter  */ 

*  Pushdown  stack  holds  two  J -bit  fields;  depth  256  */ 

*  Pushdown  stack  holds  three  fields;  depth  256  'I 

*  Pushdown  stack  holds  two  fields;  depth  256  *l 

*  Pushdown  stack  holds  one  field;  depth  256  '/ 

*  Last  value  popped  from  P] _Stackl  'I 

'  Last  value  time  popped  from  PlJStack2  'I 

*  Circular  queue  */ 

*  9  bits:    counter  pointing  to  front  of  queue  */ 

*  9  bits:   counter  pointing  to  end  of  queue  */ 

*  1  bit:  current  pixel  of  current  row  */ 

*  1  bit:  previous  pixel  of  current  row  *l 

*  /  bit:  current  pixel  of  previous  row  '/ 

*  1  bit:  previous  pixel  of  previous  row  *l 

*  2  bits:   current  brackets  *l 

*  2  bits:  previous  brackets  *l 

*  1  bit:  flag  used  in  PASS  1  V 

*  1  bit:  flag  used  in  PASS  1  V 

'flag:    indicates  down  loading  in  progress  */ 

*  flag:    indicates  Pass!  in  progress  *l 

*  flag:    indicates  Pass2  in  progress  */ 

*  flag:    enables  process  lnitialization_of_Sweep_of_Passl  '/ 
'flag:    enables  process  P2S1  Jiowjnitialization  */ 
'flag:    enables  process  P2S2_RowJnitialization  'I 
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Start_Row_Sweep_of_Passl; 

Start_Row_Sweepl_of_Pass2; 

Start_Row_Sweep2_of_Pass2; 

Next_Row_Pl; 

Next_Row_P2; 

Next_PixeLl; 

Next_Pixel_2; 

Next_Pixel_3; 

Pl_pixel_op_enable; 

P2Sl_pixel_op_enable ; 

P2S2_pixel_op_enable; 

component_number; 
current_component; 


process  Read_Pixet  'I 

process  Read_Pixel  */ 

process  ReadJPixel  */ 

process  lnitialization_of_Sweep_of  Pixel  */ 

process  P2S1  _Row_lnitialization  'I 

process  Read_Pixel  */ 

process  Read_Pixel  *l 

process  Read_Pixel  */ 

process  Pixel_Computation_of_Passl  */ 

process  Pixel_Computation_of  Sweep _of_Pass2 

process  Pixel_Computation_of_Sweep _of_Pass2 


*  Next  new  component  number  *  I 

*  component  number  of  current  scanned  pixel  *l 


'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

'  flag: 

enables 

B.2  Description  of  Primitive  Operations 

1.  Assigning  a  value  to  a  simple  variable  or  indexed  variable  has  obvious  meaning. 

2.  Four  stack  operations  are  used: 


push(stack_name,  vector); 

pop(stack_name) ; 
top(stack_name)  =  vector; 
top(stack_name)[i]  =  variable; 


/•  push  a  vector  onto  multi-field 
stack  •/ 

/*  pop  the  top  of  a  stack  */ 

/•  modify  the  top  of  a  stack  *l 

I*  modify  a  specific  field  of 
a  stack  *  I 


3.  Queue  operations  are  treated  as  indexed  accesses  to  an  underlying  array. 
B.3  Processes 

/*  Initialize  row  counters  for  Pass  I  then  start  row  sweeping  process  */ 
Process  PASSUnitialization 

if  (Passl  &  !In_Down_Load  &  !In_Passl  &  !In_Pass2)  tlien 

In_Passl  =  true; 

Row_Countl  =  511; 

Row_Count2  =  0; 

Start_Passl  =  true; 
else 

Start_passl  =  false; 
Endlf; 
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/*  Terminate  Pass  1  after  all  rows  have  been  scanned  */ 

Process  PASSl_Termination 

if  (In_Passl  &  Next_Row_Pl  &  Row_Countl=  =511)  then 

In_Passl  =  false; 
Endlf; 


/*  Initialize  column  counter  variables  and  flags  at  entrance  of  a  row  sweep 
In  Pass]  then  enable  reading  pixel  process  */ 

Process  Initialization_of_Sweep_of_Passl 

if  (Start_Passl  |  (Next_Row_Pl  &  Row.Countl  !=  511))  then 

Ncxt_Row_Pl  =  false; 

Column_Count  =  511; 

rcsct(Pl_staclcl); 

resct(Pl_stack2); 

CPIX  =  0; 

PPIX  =  0; 

touch_p  =  0; 

touch-c  =  0; 

Bracket_Countl  =  0; 

Bracket_Count2  =  0; 

Start_Row_SweepJPassl  =  true; 
else 

Start_row_Sweep_Passl  =  false; 
Endif; 


/•  At  the  end  of  each  sweep,  update  row  counters,  and  start  sweep  on  next  row  */ 
Process  End_of_Swcep_of_Passl 

If  (Ncxt_Pixcl_l  &  Column_Count=511)  then 

Row_Countl  --; 

Ro'w_Count2  --; 

Ncxt_Row_Pl  =  true; 
Endif; 


/*  Load  next  pixels  from  current  row  and  previous  row  simultaneously  into  registers  CPIX 
PPIX  respectively,  and  old  contents  of  CPIX  and  PPIX  was  backed  into  IjCPX  and  IJ'PDC 
respectively.    If  a  run  started  In  a  previous  row,  the  corresponding  bracket  marking 
is  loaded  from  bracket  table.  'I 

Process  Read_Pixel 

/*  In  sweep  1  of  pass  1 : 
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Read  pixels  from  two  rows  simultaneously  from  right  to  left. 
'! 
If  (S:art_Row_Sw(:ep_of_Passl  |  Ncxt_Pixel_l  &  Column_Count  !=  511)  then 

Ncxt_PixeI_l  =  flase; 

LCPIX  =  CPIX; 

LPPIX  =  PPIX; 

CPIX  =  image_buffer(Row_Countl,Column_Count); 

PPIX  =  imagc_buffer(Row_Count2,Column_Count); 

Column_Count  --; 

If  (PPIX  =  1  &  LPPIX  =  0)  then 

Lbrkt[left]  =  brkt[lcft]; 

Lbrkt[right]  =  brkt[right]; 

Bracket_Count2  +  +  ; 

brktfleft]  =  Lcft:_bracket_table(Row_Count2,Bracket_Count2); 

brkt[right]  =  Right_bracket_table(Row_Count2,Bracket_Count2); 
Endlf; 

If  (CPIX  =  1  &  LCPIX  =  0)  then 

Bracket_Countl  +  +  ; 

Left_bracket_table(Row_Countl,Bracket_Countl)  =  0; 

Right_bracket_table(Row_Countl,Brackct_Countl)  =  0; 
Endlf; 

Pl_pixeLop_enable  =  true; 

/*  In  sweep  1  of  pass  2: 

Read  pixels  from  two  rows  simultaneously  from  right  to  left. 
*l 
else  If  (Start_Row_SweepLof_Pass2  |  Next_PixeL2  &  Column_Count  !=  511)  then 

Next_Pixcl_2  =  false; 

LCPIX  =  CPIX; 

LPPIX  =  PPIX; 

CPIX  =  image_buffcr(Row_Countl,Column_Count); 

PPIX  =  image_buffer(Row_Count2,Column_Count); 

Column_Count  --; 

If  (PPIX  =1  &  1JPPIX  =  0)  then 

Queue_Count2  +  +  ; 
Endlf; 

if  (CPIX  =1  &  LCPIX  =  0)  then 

Queuc_Countl  +  +  ; 

Brackct_Count2  +  + ; 

brkt[left]  =  Lcft_bracket_tablc(Row_Countl,Brackct_Count2); 

brkt[right]  =  Right_bracket_tablc(Row_Countl,Bracket_Count2); 
Endlf; 

P2Sl_pixcLop_cnablc  =  true; 

/*  In  sweep  2  of  pass2: 

Read  pixels  from  two  rows  simultaneously  from  left  to  right. 

'I 
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else  If  (Start_Row_Sweep2_of_Pass2  |  Ncxt_Pixel_3  &  Column.Count  !=  0)  then 
Start_Row_Swccp2_of_Pass2  =  false; 
Ncxt_Pixel_3  =  false; 
LCPIX  =  CPIX; 

CPIX  =  image_buffer(Row_Countl,Column_Count); 
Column_Count  +  +  ; 

If  (CPIX=  1  &  1_CPIX  =  0)  then 

brkt[left]  =  Left_bracket_table(Row_Countl,Bracket_Count2); 

brkt[right]  =  Right_bracket_table(Row_Countl,Bracket_Count2); 

Bracket_Count2  --; 
Endif; 


P2S2_pixel_op_cnable  =  true; 


Endlf; 


/*  this  Is  the  main  computation  process  of  Passl.    In  different  cases  of  combinational 
conditions  appropriate  manipulation  are  carried  out  on  t\vo  stacks  and  bracket  table.  */ 

Process  PixeLComputation_of_Passl 

If  (Pl_pixel_op_enable)  then 

Pl_pLxeLop_cnable  =  false; 

case  of 

/*  Start  a  run  on  both  current  and  previous  rows  'I 
(PPIX=1  &  1_PPIX  =  0  &  CPIX=1  &&  LCPIX  =  0): 
case  of 
(brkt[rigbt]  =  0): 

push(Pl_Stackl,(l,l)); 
push(Pl_Stack2,(bracket_Countl,bracket_Countl)); 

(brkt[right]=l): 

If  (top(Pl_Stackl)[01)  then 

Left_brackct_table(Row_Countl,top(Pl_Stack2)[0])  =  1; 

Right_bracket_table(Row_Countl,Bracket_Countl)  =  1; 
else 

push(Pl_Stack2,Bracket_Countl,Bracket_Countl); 

If  (brkt[left]  =  l)  then 

top(Pl_Stackl)[0)  =  1; 

Endif; 
Endlf; 

Endcase; 

/•  Only  start  a  new  run  In  previous  row.  There  Is  no 

transition  on  current  row.  'I 
(PPIX=1  &  1_PPIX  =  0  &  !(CPIX  xor  LCPIX)): 
if  !(CPIX  &  PPIX)  then 
case  of 
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(brkt[right]  =  0): 

push(Pl_Stackl,(0,l)); 
Endcasc; 
else  case  of 

(brkt[right]  =  0): 

push(Pl_Stackl,(l,  !touch_p  |  touch_p  & 

Lbrkt[left]  =  0  &  top(Pl_Stack2)[2])); 
(brkt[nght]=l): 

If  (top(Pl_Stackl)[0]  &  (touch_p  &  !Lbrkt[left] 
&  Sl_popped[l]  I  Itouch)  then 
S2_poppcd  =  pop(Pl_Stack2); 

Rightbrackct_table(Row_Countl,S2_poppcd[l])  =  1; 
lcft_bracket_table(Row_Countl,top(Pl_Stack2)[0])  =  1; 
top(Pl_Stack2)[0)  =  S2_popped[0] 
Endif; 
Endcase; 
Endif; 

/*  Start  a  new  run  in  current  row.  There  is  no  transition 

on  previous  row.  */ 
(CPIX=  1  &  LCPIX=0  &  !(PPIX  xor  LPPIX)): 
If  !(CPIX  &  PPIX)  then 

push(Pl_Stack2,(bracket_Countl,bracket_Countl,l); 
else 

case  of 
(brkt[right]  =  0): 

If  (!touch_c)  then 

push  (bracket_Countl,Bracket_Count  1,1); 
top(Pl_Stackl)  =  (l,l); 
else 

Left_brackct_table(Row_Countl,top(Pl_Stack2)[0])  =  1; 
Right_bracket_tablc(Row_Countl,Bracket_Countl)  =  1; 
top(Pl_Stack2)[0]  =  Bracket_Countl; 
Endif; 

(brkt[right]  =  l): 

if  (touch_c  I  !touch_c  &  top(Pl_Stackl)[0])  then 

Right_brackct_tablc(Row_Countl,Bracket_Countl')  =  1; 

Lcft_bracket_tablc(Row_Countl,top(Pl_Stack2)[0])  =   1; 

top(Pl_Stack2)[0]  =  Bracket_Countl; 
else  If  (!touch_c  &  ltop(Pl_Stackl)[0])  then 

push(Pl_Stack2,Bracket_Countl,Bracket_Countl); 

if  (brkt[left]=l)  then 

top(Pl_Stackl)[0]  =  1; 

Endif; 
Endif; 

Endcase; 
Endif; 

/*  Reached  the  end  of  a  current  run  in  both  previous  row 

and  current  row.  */ 
(PPIX=0  &  1_PPIX=1  &  CPIX  =  0  &  LCPIX=1): 


B.6 


Appendix  B  CCB 


touch_c  =  0; 

touch_p  =  0; 
case  of 
(brkt[left]  =  0): 

if  (top(Pl_Stackl)[l])  pop(Pl_Stack2); 

Sl_popped  =  pop(Pl_Stackl); 
Endcase; 

/*  Reached  the  end  of  a  run  In  previous  row.  But  there  is 

no  transition  on  current  row.  */ 
(PPIX  =  0  &  1_PPIX=1  &  !(CPIX  xor  LCPIX)): 
If  (!CPIX)  then 
case  of 
(brkt[lcft]  =  0  &  brkt[right]  =  0): 

If  (touch_c  &  top(Pl_Stackl)[l])  then 

pop(Pl_Stack2); 
Endlf; 
Sl_popped  =  pop(Pl_Stackl); 

(brkt[lcft]  =  0  &  brkt[nghtl  =  l): 

If  (top(Pl_Stackl)[0]  &  top(Pl_Stackl)[l]  | 
touch_c  &  top(Pl_Stack2)[2])  then 
pop(Pl_Stack2); 
Endlf; 

Sl_poppcd  =  pop(Pl_Stackl); 
Endcase; 
elie 

touch_p  =  1; 
case  of 
(brkt[left]  =  0): 

top(Pl_Stack2)[2I  =  top(Pl_Stackl)[l]; 
Sl_poppcd  =  pop(Pl_Stackl); 
Endcase; 

Endlf; 

touch_c  =  0;  , 

/*  Reached  the  end  of  a  run  in  current  row.  But  there  is 

no  transition  on  previous  row.  */ 
(CPIX  =  0  &  LCPIX=1  &  !(PPIX  xor  LPPIX)): 
If  (IPPIX)  then 

If  (!touch_p  I  touch_p  &  brkt[left]  =  0  &  top(Pl_Suck2)[2])then 

pop(Pl_Stack2); 
Endlf; 
else 

touch_c  =  1; 
Endlf; 
touch_p  =  0; 

/•  Start  a  new  run  in  previous  row  and  reach  the  end  of  the 

current  run  in  current  row  simultaneously .  */ 
(PPIX=1  &  LPPIX  =0  &  CPIX=0  &  LCPIX  =1): 

If  (!touch_p  I  touch_p  &  brkt[lcft]  =  0  &  top(PLStack2)[2])then 
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pop(Pl_Stack2); 

Endlf; 

If  (brkt[right]  =  0)  then 

push(Pl_Stackl,(0,l)); 
Endlf; 

touch_p  =  0; 
touch_c  =  0; 

/'  Start  a  new  run  In  current  row  and  reach  the  end  of  the 

current  run  In  previous  row  simultaneously .  *l 
(CPIX=1  &  LCPIX  =  0  &  PPIX  =  0  &  1_PPIX=1): 

If  (brkt[left]  =  0  &  touch.c  &  top(Pl_Stackl)[l])  then 

pop(Pl_Stack2); 
Endlf; 

push(Pl_Stack2,Bracket_Countl,Bracket_Countl,l); 
touch_p  =  0; 
touch_c  =  0; 

Endcase; 

Next_Pixel_l  =  true; 
Endlf; 


/•  Initialize  row  counters,  queue  counters  and  set  next  component  number 
to  1 ,  then  start  first  sweep  of  row  1  */ 

Process  PASS2_Initialization 

if  (Pass2  &.  !In_Down_Load  &  !In_Passl  &  !In_Pass2)  then 

In_Pass2  =  true; 

Row_Countl  =  0; 

Row_Count2  =  511; 

Queue_Countl  =  0; 

Queue_Count2  =  0; 

component_number  =  1; 

Start_Jlow_Sweep2_of_Pass2  =  true; 
else 

Start_Row_Swccp2_of_Pass2  =  false; 
Endif; 


/*  Terminate  pass  2  after  all  rows  have  been  scanned  */ 

Proceis  PASS2_Tennination 

if  (In_Pass2  &  Next_Row_P2  &  Row_Countl=0)  then 

Iii_Pass2  =  false; 
Endif; 
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/•  Initialize  column  counter,  stacks  and  other  variables  for  the  first 

sweep  of  a  row  In  pass  2;  then  enable  process  Read_Plxel  to  load  next 
pixels  from  current  row  and  previous  row  */ 

Process  P2Sl_RowJnitialization 

if  (Start_Pass2  |  (Ncxt_Row_P2  &  Row.Countl  !=  511))  then 

Ncxt_Row_P2  =  false; 

Column_Count  =  511; 

reset(P2_stackl); 

rcset(P2_stack2); 

CPIX  =  0; 

PPIX  =  0; 

Bracket_Count2  =  0; 

Start_Row_Swccpl_of_Pass2  =  true; 
else 

Start_Row_Sweepl_of_Pass2  =  false; 
Endif; 


/•  After  the  first  sweep,  from  right  to  left,  reaches  the  end  of  a  row, 
second  sweep ,  from  left  to  right,  is  enabled  '! 

Process  End_of_Sweepl_of_Pass2 

If  (Next_PixeL2  &  Column_Count  =  511)  then 

Sweep2_of_Pass2_cnablc  =  true; 
Endlf; 


/'  This  Is  main  computation  process  In  first  sweep  of  a  row 
In  pass  2  'I 

Process  Pixel_Computation_of_Swcepl_of_Pass2 

If  (PlSl_pixcl_op_enable)  then 

PlSl_pixeLop_cnablc  =  false; 

case  of 

/•  Start  a  new  run  on  both  previous  row  and  current  row 

simultaneously.  'I 
(PPIX=1  &  1_PPIX=0  &  CPIX=1  &&  LCPIX=0): 
case  of 
(brkt[right]  =  0): 

push(P2_Stackl,(Queuc(Queue_Count2),l)); 
(brkt[right]=l): 

if  ltop(P2_Stackl)[l)  then 

push  (P2_Stack  1 ,  (Qucue(Queue_Count2) ,  1) ) ; 
Endlf; 
Endcase; 
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/*  Start  a  new  run  in  previous  row .  But  there  is  no 
transition  on  current  row.  */ 
(PPIX=  1  &  LPPIX=0  &  !(CPIX  xor  LCPIX)): 

If  (CPIX  &  PPIX  &  !top(P2_Stackl)[l])  then 

top(P2_Stackl)  =  (Qucuc(Qucue_Count2),l); 
Endlf; 

/*  Start  a  new  run  in  current  row.  But  there  is  no 

transition  on  previous  row .  */ 
(CPIX=1  &  LCPIX  =  0  &  !(PPIXxorl_PPIX)): 
case  of 
(brkt[nght]  =  0): 

If  (CPIX  &  PPIX)  then 

push(P2_Stack(0,0)); 
else 

push(P2_Stackl,(Qucue(Queuc_Count2),l)); 
Endlf; 
(brkt[right]  =  l): 

If  (CPIX  &  PPIX  &  !top(P2_Stackl)[l])  then 

top(P2_Satackl)[l]  =1; 
Endlf; 
Endcasc; 

/*  Reach  the  end  of  the  current  run  on  both  previous 

row  and  current  row  simultaneously .  '/ 
(PPIX=0  &  LPPIX=1  &  CPIX  =  0  &  LCPIX=1): 
If  (brkt[left]  =  0)  then 

push(P2_Stack2,pop(P2_Stackl)[0]); 
Endlf; 

/*  Reach  the  end  of  current  run  in  previous  row.  But 

there  Is  no  transition  on  current  row.  */ 
(PPIX=0  &  1JPPIX=1  &  !(CPIX  xor  LCPIX)): 
NOP; 

/•  Reach  the  end  of  current  run  In  current  row  and 
there  Is  no  transition  or  start  a  new  run  In 
previous  row.  */ 
(CPIX  =  0  &  LCPIX  =  1  &  !(PPIX  xor  LPPIX)) : 
(PPIX=1  &  LPPIX=0  &  CPIX=0  &  LCPIX=1): 
case  of 
(brkt[lcft]  =  0): 

if  (top(P2_Stackl)[l])  then 

push(P2_Stack2,pop(P2_Stackl)[01); 


else 


Endlf; 
Endcase; 


push(P2_Stack2,componcnt_numbcr); 
pop(P2_Stackl); 
componcnt_number  +  +  ; 


/*  Start  a  new  run  In  current  row  and  reach  the  end 
of  current  run  in  previous  row.  */ 
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(CPIX=1  &  1_CPIX  =  0  &  PPIX  =  0  &  1_PPIX=1): 
case  of 
(brkt[right]  =  0): 

push(P2_Satackl,(0,0)); 
Endcase; 

Endcasc; 

Ncxt_Pixel_2  =  true; 
Endlf; 


/•  Initialize  column  counter,  stack,  and  other  variables  for  second 
sweep  of  a  row  in  pass  2;  then  enable  process  Read_Plxel  */ 

Process  P2S2_RowJnitiaJization 

If  (Swecp2_of_Pass2_cnable)  then 

Swcep2_of_Pass2_cnable  =  false; 

Column_Count  =  0; 

rcsct(P2_stackl); 

CPIX  =  0; 

current_component  =  0; 

Qucue_Count2  =  Queue_Countl; 

Start_Row_Swcep2_of_Pass2  =  true; 
cndif; 


/*  After  the  second  sweep  of  a  row  is  finished,  row-counters  are  incremented 
by  one  for  next  row.  */ 

Process  End_of_Swcep2_of_Pass2 

If  (Next_Pixel_3  &  Column_Count=0)  then 

Row_Countl  +  +  ; 

Row_Count2  +  +  ; 

Ncxt_Row_P2  =  true; 
Endlf; 


/•  This  Is  main  computation  process  In  second  sweep  of  a  row  In  Pass  2  */ 

Process  Pixel_Computatioii_of_Swccp2_of_Pass2 

if  (PlS2_pixel_op_enable)  then 

PlS2_pixeLop_cnablc  =  false; 

case  of 

/*  Start  a  new  run  In  current  row.  */ 
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(CPIX=1  &  1_CPIX=0): 
case  of 

(brkt[left]  =  0  &  brkt[right]  =  1): 

current_component  =  pop(P2_Stack2); 

push(P2_Stackl,currcnt_componcnt); 
(brkt[left]  =  l  &  brkt[right]=  1): 

currcnt_component  =  top(P2_Stackl); 
(brkt[lcft]  =  l  &  brkt[right]  =  0): 

currcnt_component  =  pop(P2_Stackl); 
(brkt[lcft]  =  0  &  brkt[right]  =  0): 

current_component  =  pop(P2_Stack2); 
Endcase; 

Queuc(Queue_Count2)  =  currcnt_componcnt; 
Queue_Count2  +  +  ; 

/*  Reached  the  end  of  current  run  in  current  row.  */ 
(CPIX=0  &  LCPIX=1): 

current_component  =  0; 

Endcase; 

output(current_component); 

Next_Pixel_3  =  true; 
Endif; 
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